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ABSTRACT: Thermodynamic data for nitrogen are presented in 
Mollier diagram Yorm covering the temperature range of 30°K 

ee ot to 15,000°K and the density range of 10°” to 20° amagat. Data 
ae for the selid and liquid regions are included. ‘The diagram 

: contains all available data of interzst to the worker in fluid 
mechanics and aerodynamics. Speed of sound data are prssented 
on separate charts. 
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Mollier Diagram for Nitrogen 


This report brings together in one large scale diagram most 
of the available data for nitrogen. The temperature range 
covered is 30°K to 15,000°K while densities from 197” to 10° 
are represented. This diagram is especially recomended fcr 
aerodynamicists and those engaged in aerodynamic facility 
design. 


The authors acknowledge the assistance of Mr. John R. Bott 
who did the lettering of the diagram sheets. 


E. F. SCHREI' ER 
Captain, USF 
Commander 
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Hot ca ANE 


INTRODUCTION 


The dezand for high Mach number low altitude flight simulation 
in wind tunnels has increased in recent years. Wind tunnel supply 
temperatures and pressures have increased tremendousiy. Tunp<sis 
now being planned require pressuresto 5000 atmospheres and teupsra- 
tures to 2800°K. At conditions such as these there are tide cevis- 
tions fros the ideal gas and such things as wariable specific heat 
and interszolecular forces must be considered. 


The need for adequate thermodynamic property data under these 
conditions is evident. The past ten years have seen an outpouring 
of works on the thermcdynamic properties of air and these properties 
mzy now be considered to be weli documented. Hcwever, the study of 
the properties of nitrogen has only recently resulted in adequate 
tables for high temperatures. 


Thermodynamic dizgrams have several advantages over tabuiated 
tables. One of the principal advantages is convenience. The dia- 
gram allows trends to be seen at a glance and presents masses cf 
data in a few sheets. The required information on mcst common 
thermodynamic processes can be obtained from the diagram without 
tedious interpolation in tables. In most cases the accuracy of 
the diagram is sufficient for engineering use. 


The making of thermodynamic diagrams for gases has been a 
popular activity among thermodynamicists and chemists for the past 
60 to 70 vears. Din (ref. (1)) in his three volume work on gases 
for industrial use gives an excellent historical and critical 
review of such diagrams in general (ref. (la)) and for particular 
gases as he discusses them individually. Generally, industrial and 
chemical users have preferred the temperature-entropy diagrasz, for 
example references (1) and (2), while aerodynamicists and power plant 
engineers used the Mollier or enthalpy-entropy diagram (refs. (3)- 
(7)). The many isenthalpic and isentropic processes encountered in 
fluid dynamics make the Mollier diagram popuiar in this field. 


Yost of the existing diagrams for nitrogen are limited in 
their coverage. Almost all do not include the high density-high 
temperature region that is now of growing interest. A recent report 
by Grabau and Brakinsky (ref. (%)) gives calculated properties of 
nitrogen for densities to 1000 amagat. This filled in a region for 
which no adequate data had been available. 
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The large smount of recent data that is now availabie and the 
limited sccpe of other diagraszs made it advisable to cozpile a new 
diagram which would include as much of the available data as possi- 
ble. It also afforded the opportunity to compare and study the 
waricus sources of thermodynamic properties. The new diagra= is, 
presented in this report. It covers a range froz 30°K to 15,000° 4 
and a density range from log p/p,= -7.0 to 3.0; data for the solid, 
liquid, and wapor regions have been inciuded. 


Tbe sources, the compilation, and the use of the diezgranm will 
be discussed. Appendices containing information helpful ir the 
use of the diagram have been inciuded. 


SODRCES FOR THE DIAGRAM 


A great many sources were used in the preparation of the diagraz. 
Usually each source covered only a limited region. In some cases 
units and reference points had to be converted. Figure 1 gives a 
skeleton chart shoring the scope cf the diagra=s 2nd reference nuzbers 
for each region. In some areas the regions overlap and several 
sources were used. 


The bulk of the data from 100 to 1500°K and a density range 
from log ¢/p,= -7.0 to 2.8 was taken from Little and Neel (ref. (9)). 
Their tables" incorporate the work of Din (ref. (1Ic))} and Hilsenrath, 
et al (ref. (10)) for pressures down to 10°* atsospheres. They used 
a perfect gas extrapolation for pressures beiow 10°? atmospheres. 

At the high densities, the tabuiated pressures range fron 16* ateos- 
pheres at 300°K to 10° atmospheres at 1500°K. 


The tables of Din are based on the best available experizental 
data. Din believes that his worst points are no more than 1.5 per- 
cent in error while most of the tables are good to within 0.5 percent. 
Din made an extensive comparison of the work of many experimenters 
and from these chose the best data for calculating his tables. 


Below 100°K the ideal gas thermodynamic functions for zolecular 
nitrogen of Hilsenrath, et al (ref. (1)) were used to detersine 
enthalpy levels at zero pressure and entropy values at stazospheric 
pressure. Temperature lines were found by assuming that (@H/3p) = 0. 
The values of entropy at low pressures were found from the ied 
perfect gas relationship. 


(S/R), = (G/2). - inp 
p= 1 atm 
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The iine co.mecting the critical point and the triple point 
was taken from Din (ref. (ic)). Some pressure and density lines 
near this region were taken fron reference (2). This region 
defines the bounds between the saturated vapor and the caturated 
liquid. it is equivalent to the zero soisture line on 2 steaz 
Mollier diagram. From the triple point to lower pressures the line 
of equilibrium between the sclid and vapor phases (zero vapor line) 
was Gefined by the expression from Din (ref. (ic)) converted fro 
mm Eg to atzospheres: 


= (255.821) 
= Zz == 
iog Pp 3.61513 = S-c05) 


where T is in °Z. 


A discrepancy exists between the ideal gas data of reference (10) 
and the experimental triple point obteined fro= Din (ref. (ic)) . The 
difficulty has not been resolved and caution should be used in inter- 
preting the diagram near the zero vapor boundzry beicw 100°K. Away 
from the bouncary the ideal gas data should be adequate. 


Data for the smelting soiid and the freezing iiquid were taken 
from Dis (ref. (ic)) and frow reference (2). Extrapolations vere 
made to extend these data to 1000 amagat Gor p/p,= 3.0). Properties 
at pressures above the critical point were taken from in (ref. (2)) 
and from Little and Neel (ref. (9)}. The trend of the data in this 
ragion indicates that a point is reached about 1600 azagat shere the 
solid and the liquid regions coincide. 


Properties for the high temperature lor censity region were 
taken from the tables of Hiisenrath and Kiein (ref. {11)). These 
tables include the effects of dissociation, ionization, and inter- 
molecular forces. intermolecular forces are included to the seccznd 
virial correction. Hiisenrath and EKiein use the weii-xnows procedure 
of listing aii of the isportant chemical teactions takizg piace, 
determining the ancunt of reactants prercent at each condition and 
calculating the therzodyramic proverties of the zwixture. These 
sroperties are then corrected Zor the interzoleczvisr forces. ‘The 
tables have been widely used and zre accepted 2s reljahie. 


Grabau and Brahinasxy {ref. (8)) recentiy comilied 2 set of 
tables for high density nitrogen {to i000 azagat) which includes 
third virial corrections. These tables extend to 5000°Z and hare 
been designed to biend well with the tables of Littie and Heei 
(ref. (8)) 2t low tezperateres. At the higher temceratures the 
property values of Hilsenrsth and Klein (ref. (11)), and Ssith 
(ref. (312)) at log p/p,= 0 verte used 2s points of departure. While 
Grabau makes 20 claiss for accuracy, he feeis that the tables sre 
internally consistent and are adequate for engineering use. 


‘a 
wh 
1 


bey PAPER ETE NNN THGS 70 99h UEPU LENA EP SPALL DCE PUREE bags cece tee 


gute co ga Up onyas ovtRbeyenae pAb a tLe TreetooeffAPbeeaetez tag eat + 


OT 


Ton 


punter 


E 
| 
el 
= 


ean 
“i 
i 


ELTVANCOU AVENE REE PEUELGERU LARGE sPovabianELeroeReROCepeOn EAL HARSHA 


aE hE te alec pen 


iT 


PT Ld) 


‘fen + 


My ete ae 


EGLTR 67-32 


COMPILING THE NITROGEE MOLLIER DIAGESY 


The convenient forust of Littie (ref. (3)) and Eowhrey, Littie, 
and Seeley (ref. (4)) was chosen for the diagraz. The disensicaiess 
forzs, such as H/RT (enthalpy, S/2 (entropy), log p/p, {density in 
anarat units}, and P (ateospheres) used in the Eatiozai Bureau of 
Standards (285) tables (refs. @0) acd {11)) and those based on the 
EBS work jrefs. GG), (4, (6), (8), 2nd (8)) were used ehenever 
possible. The enthalpy faenction is shoen as H/2 in units of "E. 
Conversicn factors for various dizessicn systees ar : 


a 
i 
b 


log p/p, is given in Appendix C. 


Sections of the dizgra= are plotted on sexilog paper. Each 
sheet covers one order of sasnitede of enthaipy and seren units of 
entropy. 4n index to the sections is given in figure 2. This #311 
aid ini tins the desired sections. The diagram is cat off at 2 
v2inue of S/2 = 49. The area of serious interest to workers in fiuic 
Gynazics lies in the region there ralces of entropy are low. A loot 
at figure 2 shows that the diagras is beginning to taper to a ioxng 
point bounded by FT = 15,600°K and log 5/p,= -7.9. Thus only a small 
avount of useful informations is lost by stopping the diagraz 2t 
S37 = 48, 


Whenever the work of Din (refs. (ic) asd (2)> eas used directi;, 
the tabclated values were charged to tie dimezsionless vaiues of BS. 
The different reference points for zero entropy a=d enthaipy were aiso 
adjusted. The equations used are giren in Appendix DB. 


The interpolation necessary to fill the low cessity regic=a 
between 1500 and 2000°K was dome graphically. Slight differences 
between tables sere smocthed by fairizcg 2 carve through sererai 
points on the curve. 


GSiSG TSE ZGLLIER DIsAGGAR 


This section is éesisgned to heip the user of the diszTaz 2nd 
may be considered as independent of the remainder of the text. 


The Mollier dizsraz contains waiues for the pressure, density, 
tenperatore, enthalpy, and entropy of aitrogen. Temperature ("EK is 
represented by heavy sciid lines and pressure (in ateospheres) by 
thin selid lises. The density (ice s/z,, pp in amagat units} is 2 
dashed line. Enthaipy is in units of sc the entrosr is disensics- 
jess. Pressure, temperature, and density lines are spaced at cicse 
The syxbols used are given in Appendix 4 and coaversion factors to 
obtain other units in varices dimension syste=s are in Apsencix B. 
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The compressibility factor, Z, may be found from the expression, 
% = P/pRT. This is especially easy as R is 3.66i x 1073 (1/(273.16°K)) 
when F is in atmospheres, p_in amagats, and Tin °K. If T is in °R, 
R is 2.175 x 107*°(1/(459.69°R)). These values of R are independent 
of the type of gas so long as atmosphere and amaget units are used. 


The index sheet, figure 2, will help to locate the proper sheets 
within the diagram. The sheets are numbered in columns so that an 
isentropic process may te followed more easilvr. 


SOURCES FOR THE ACOUSTIC VELOCITY CHARTS 


The speed of sound data were taken from Lewis and Neel (refs. (13) 
and (14)). They calculated specific heats and acoustic velocities 
from the tables of Little and Neel (ref. (9)), Hilsenrath and Klein 
(ref. (11)), and Smith (ref. (12)). The tables were plotted at 
constant density as a funct‘.on of temperature to check the matching 
of the various sources. The acoustic velocity blended smoothly from 
source to source except at the very high densities. At densities of 
log p/p, = 2.4 and 2.6 discontinuities exist between acoustic veloci- 
ties calculated from Little and Neel and those calculated from Smith. 
The discontinuities are modest at log p/p,= 2.4 but are very large at 
log p/p,= 2.6. The c: 2xts probably should not be used for densities 
greater than log p/p, 2.4. The acoustic velocities are given in 


terms of a/e, where a, is tke acoustic velocity at standard conditions. 


Vaiues of a, are given in Appendix B. The curves are given at con- 

stant temperature as a functien of log p/p,. The three charts cover 
the low temperature (T < 2000°K)-low density region (log p/p, <9), 

the high temperature region (T > 2000°K), and the low temperature 

(T < 2000°K)-high density region (log p/p, > ~1.0). 
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ae 


a acoustic -velocity 
: a, acoustic velocity, S.T.P. 
°c temperature in degrees Centrigrade 
Cop. critical point 
°F temperature in degrees Fahrenheit 
H enthalpy 
°K temperature in degrees Kelvin 
Pp pressure (atm) 
De reference pressure 
R gas constart 
R temperature in degrees Rankine 
Ss entropy 
T temperature 
t.p. triple point 
: PA compressibility factor 
P density 
Py reference density 
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APPENDIX B 


STANDARD PROPERTIES AND CONVERSION FACTORS 


1. Standard Properties 
mclecuiar weight 


density of the gas at S.T.P. 
(iuternational atm. g = 980.665) 


molecular volume at N.T.P. 
boiling point cf liquid at 1 atm. 


triple point: temperature 
pressure 


critical constants: temperature 
pressure 
volume 
density 


vf THUS THAR TENT Er TEE HHH 


Shi MB emmy Ag 


VAD Tp) 


2. Convernsion Factors 


juatinny 


28.016 


G.0012525 grams/cm® 
22,403.5 cm® 
77.35°K 


63.15°K 
94.01 mm Hg 


126 .2°K 

33.5 atm. 
90.094 cm*/mole 
248.668 amagat 


Standard Atmospheric Pressure Reference De-sity 


Pog ftingyiee 


P, = 760 mm Hg Pe = 1.25046 x 107° gram/cm® 
29.921 in He 4.46338 x 107° moles/cm® 
10332 kg/m? 1.25050 gram/liter 
- 14.696 psia 4.51760 x 1075 1b/in® 
2116 psfa 7.80462 x 10°* lb/ft? 
2.42435 x 107° slugs/ft® 


Temperature Acoustic Velocity 
T (°K) x1.8.=T (°R) a, = 1105.66 ft/sec 
T ((K) = T Co + 273.16 = 0.337006 mm/psec 
T CR) = T (FP) + 459.69 = 0.337006 km/sec 
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With Units 
To Change To of Multiply By 
B/R (°K) H/RT none 1/T (°K) 
H/R °R 1.8 
H £t? /sec? 3196.66 
BTU/1b 0.12768 
cal/gm 0.0709305 
S/R s ft? /sec?-°K 3196.66 
ft? /sec?-°R 1775.92 
BTU/1b-°R 0.07093 
BTU/1b-°K 0.12678 
” eul/gm-°K 0.0709305 
p/p, ? gm/cn? 1.25046 x 107° 


Lb/in® 4.51760 x 107° 
lb/ft? 7.80641 x 107? 
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| APPENDIX C 
" DENSITIES CORRESPONDING TO THE DIAGRAMMED ARGUMENTS 
zs . Log p/p, p/p, (amagats) 
z awe ere ii i eat, Sach ode 
z 3.0 1000.000 
= 2.9 794.333 
| 2.8 630.957 
Ed 2.6 398.107 
4 2.4 251.189 
s 2.2 158.489 
a 2.0 100.000 
5 1.8 63.0957 
= 1.6 39.8107 
= 1.4 25.1189 
; 1.2 15.8489 
¥ 1.0 10.0000 
z 0.8 6.30597 
z 0.6 3.98107 
: 0.4 2.51189 
‘ 0.2 1.584489 7 
3 0.0 1.000000 4 
t -0.2 0.630957 
5 -0.4 0.398107 
A -0.6 0.251189 
: -0.8 0.158489 
3 -1.0 0.100000 
5 -1.2 0.0630957 
4 ~1.4 0.0398107 
i -1.6 0.0251189 
-2.0 0.0100000 
: etc. etc. 
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APPENDIX D 


CONVERSION FACTORS FOR F. DIN (REF. 


Din's units are: 


Entropy - joules/(moie-°K) 


Volume - cubic centimeters/mole 


Enthalpy - joules/mole 


p/p, = 2.24035 x 10*/(tabulated value) 


S/R = 1.2173 x 107: x tabulated value 


H/R 


= 


(ic}} 


= (1.2173 x 107? x tabulated value) - 848.0 
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FIG. 23 ACOUSTIC VELOCITY FOR LOW TEMPERATURE NITROGEN 
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